Time reversal Aharonov-Casher effect in mesoscopic rings with Rashba spin-orbital 
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The time reversal Aharonov-Casher (AC) interference effect in the mesoscopic ring structures, 
based on the experiment in Phys. Rev. Lett. 97, 196803 (2006), is studied theoretically. The 
transmission curves are calculated from the scattering matrix formalism, and the time reversal 
AC interference frequency is singled out from the Fourier spectra in numerical simulations. This 
frequency is in good agreement with analytical result. It is also shown that in the absent of magnetic 
field, the Altshuler-Aronov-Spivak type (time reversal) AC interference retains under the influence 
of strong disorder, while the Aharonov-Bohm type AC interference is suppressed. 
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The advancement of spintronics and spintronic ma- 
terials has provided opportunities to study and utilize 
new electronic devices based on the electron spin de- 
grees of freedom The novel spin related properties 
can be detected by tuning the spin-orbital interaction 
(SOI). Phenomena related to Rashba effects in, e.g., 
InAs and Ini-xGa^As based two-dimensional electron 
gas (2DEG) systems have been observed The value 
of the Rashba parameter in these semiconductor systems 
can be as high as 4 x 10~ n eVm. The experimental 
results @, i] have shown that in InAs and Ini^Ga^As 
based 2DEG systems, the Rashba effect is responsible 
for the spontaneous spin splitting. Many advanced elec- 
tronic devices have been proposed to manipulate electron 
spin, such as spin transistor 0, waveguides [f| Q, spin 
filters spin interferometers!^, Q , etc. In the spin in- 
terference devicefl, [t| the resistance is adjusted via 
the gate voltage using the spin interference phenomenon 
of the electronic wave functions in a mesoscopic ring. In 
the presence of Rashba SOl[ll|, the wave function ac- 
quires the Aharonov-Casher (AC) phase when an elec- 
tron moves along the ringjl^, [Hj]. If the quantum coher- 
ence length is larger than or comparable with the elec- 
tron mean free path, the resistance amplitude will depend 
on the AC phases of electron wave functions in the two 
arms of the ring and display interference pattern. Mean- 
while, the modulation of interference with the aid of AC 
phase could be further tuned through an applied gate 
volt age [3|. 

The experimental demonstrations of the spin interfer- 
ence in the AC rings have been reported recently 14 , TBI . 
As expected, the resistance of the ring oscillates with the 
gate voltage and an external magnetic field. However, 
their experimental results do not directly establish the 
time reversal AC frequency. In this study, we numeri- 
cally investigate the transport properties of multi meso- 
scopic rings with magnetic field and Rashba SOI in the 
framework of Landauer-Biittiker formalism. The inter- 



ference modes are analyzed by means of the fast Fourier 
transformation (FFT) algorithm. As a subsidiary in- 
terference, the magnetic flux leads to the conventional 
Aharonov-Bohm (AB) [l|| and Altshuler-Aronov-Spivak 
(AAS) [13, El effects. Subtracting the contribution of 
magnetic flux and determining the pathes of electrons 
from the frequency spectra, AB and AAS type oscilla- 
tions associated with SOI in rings are confirmed. In the 
presence of dephasing, AAS type (time reversal) AC in- 
terference remains, while the AB type oscillation is sup- 
pressed. 
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FIG. 1: The schematic graph of the M-N ring-array structure. 
The dot lines represent the omitted rings. 6 is the local polar 
coordinates. 

The M-N ring-array structure, as shown in Fig. [1] 
mimics the experimental configurations [lij ]. For sim- 
plicity, each arm of the ring is assumed to be a one- 
dimensional (ID) quantum wire and connected with each 
other at crossing points. The 2M leads are connected 
to the M-N ring-array as M ID-channels. The exter- 
nal magnetic and electric fields are perpendicular to the 
plane of the rings. The system can be solved by dividing 
it into several sections, i.e., the semi-infinite left M- leads, 
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the region of M-N rings array in the presence of SOI, and 
the semi-infinite right M-leads. With the continuity con- 
dition on the boundaries between sections, the scattering 
matrix formalism (l9j can be established. 

Because of the absence of the SOI in the leads, 
the wavefunction in the left (right) i-th channel can 
be written in the form of ipiL = a,iL, + e lkx \ j> 
+a lL ^e ikx \ |> +b lL ,+e~ tkx \ T> +b lL ^e^ kx \ J> and 
i>iR = a iR ^ + e lkx \ f> +a iR -e lkx \ j> for an inci- 
dent electron with momentum hk. The coefficients sat- 
isfy the condition Y^i„=± \ b iL.a\ 2 + J2ia=± l a »-R,o-| 2 = 

J2ia=± \ a iL,a\ = 1. 

In the region of M-N rings array, the electronic wave- 
function \1/ in the rings can be expanded in terms of 
the complete set of eigenfunctions ip na (9) as (9; k) — 
J2na c naipn<T (9; k) with the spin index a for the ring (i, 
j). For a single ring, the Hamiltonian in the local polar 
coordinates reads [2Cj| 
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(Here m* is the electron effective mass and we take 
TO* = 0.05m, r is the radius of a single ring, a 
is the Rashba SOI strength, and <I> is the mag- 
netic flux in unit of the flux quanta <i>o = h/e. 
We ignore the Zeeman terms because of the weak 
magnetic field, and concentrate on the effect gener- 
ated by SOI.) The eigenfunctions of the system are 



i>n,- (9) 

and the 
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corresponding eigenvalues are e n> ± ■ 
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2mra/h 2 , 

fl± = —1/2 (1 ± cos/3), and T denotes the transposition 
of vectors. 

The linear equations among coefficients c na can be 
determined by the Griffith conditions 2l| at each joint 
point. After obtaining the sets of wavefunctions tpiL, ipm 
and (9; k), the scattering matrices can be derived ac- 
cording to the multi-mode scattering-matrix procedure 
given in Ref. [19]. Therefore, the total transmission 
and reflection coefficients of ring-array system can be ob- 
tained. 

We have calculated the transmission coefficients for 
several ring-array systems under different magnetic flux 
<I> and Rashba coefficient a. The phenomena are similar 
in the systems with different number of rings and radii, 
and we take the 3-3 ring- array with r = l/_tm as a special 
case to display our results. In Fig. [D[a), we present the 
transmission coefficient versus magnetic flux (T vs $) of 
the 3-3 ring-array with the electronic incident wavevec- 
tor k = 1.009651 x 10 8 m -1 and the Rashba coefficient 



a = 2peVm. The transmission coefficient oscillates with 
the magnetic flux periodically. Although the shapes of 
the curves may look quite different with different k and 
a, the periods of the curves are the same. In fact, as 
we can see in the right panel, the FFT spectra of the 
curve reveals different interference modes (namely AB, 
AAS, etc.). For various systems with different k and a, 
the FFT spectra of the T vs $ curves include the same 
modes, but the amplitude of these modes will be differ- 
ent. 

In the absence of magnetic flux, transmission T oscil- 
lates with tuning the SOI strength [Fig. Hfb)] . The oscil- 
lation is no longer strictly periodic at large incident elec- 
tron momentum, because of the complicated back and 
forth multi-scattering between channels and rings. How- 
ever, the Fourier spectra still reveal the same main fre- 
quency at a -1 = 0.646 (peVm) , which corresponds to 
the AB type oscillation. Besides this mode, other modes 
at a -1 = 1.312 (peVm)" 1 , ... (AAS type etc.), are not 
very clear here. While in the case for a single AC ring, 
these modes will be clearly shown in the FFT spectra. 

The frequencies subsumed in the Fourier spectra on 
the right panels of Fig. [DJb) can be understood ana- 
lytically from a single ring structure. The AC phase 
difference of opposite spins traveling in full counter- 
clockwise (CCW) and clockwise (CW) circle are given 

by M A^£. 
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In the limit of 



large a, the two phase differences take the same form as 
in Ref. [a. flit] . Therefore, the AAS oscillation amplitude 
can be expressed as 22 1 
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From this relation the AAS type (time reversal) AC 
frequency is approximately at a~^ AS = 2rm*/h 2 — 
1.31peVm, while the AB type oscillation frequency is 
a~A B — rm*/H 2 = 0.65peVm, as indicated in the Fourier 
spectrum. 

In order to emphasize on the AAS effect in the ring 
structure and bring light to the objective law govern- 
ing the spin interference, we introduce disorder induced 
dephasing effect, so that the AB type oscillation is sup- 
pressed. For simplicity, a single ring with two leads is 
considered, which is enough to manifest the main princi- 
ple of AAS effect in a ring structure. The conductance is 
numerically calculated with a tight-binding model. The 
tight-binding Hamiltonian for a ring in Eq. (1) is given 



3 



(a) 





2 4 
o / [peVm] 



AB AAS 




1 2, 

a I [peVm] 



FIG. 2: (a) The curves for T vs $ and the corresponding 
Fourier spectrum for 3-3 ring-array with k — 1.009651 x 
10 8 m _1 and a — 2peVm; (b) The curves for T vs a and 
the corresponding Fourier spectrum for 3-3 ring-array with 
k = 1.009651 x lO 8 ™" 1 and $ = 0. 



oscillation amplitude remains. Similar results have been 
reached in a ID square loop system 
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To see the suppression of AB type oscillation in the 
presence of disorder, the conductance is calculated for a 
single ID ring with different Rashba SOI strengths. Em- 
ploying FFT and inverse FFT techniques, we first cal- 
culate the flux dependent conductance and then extract 
the AAS oscillation amplitudes at zero magnetic field for 
a series of different Ao- Afterwards we retrieve the in- 
formation about AAS oscillation amplitudes vs Ao (same 
steps as in Ref. [15|). The purpose of this approach is to 
single out the time reversal AC phase differences. Under 
the influence of magnetic field and Rashba SOI, two paths 
will accumulate both AB and AC phases before interfer- 
ence. Since AAS oscillation is caused by the interference 
of two time reversal paths (one CCW circle and another 
CW circle), it is feasible to find out AC phase difference 
between these two time reversal paths, if we extract AAS 
oscillation and set magnetic field to be zero. 
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where e; is the lattice on site energy, the operator 
c Jo-(q,(t) creates (annihilates) an electron with spin a at 
site I, N is the total number of lattices in the ring, and 
ipi = 7r(2Z — 1)/N as defined in Ref. 23]. The hopping 
coefficients in the second and third summations can be 
expressed as t (and A) = to (and Ao) exp(ie$/ftA r ) and 
to = h 2 /2m*a 2 is the nearest neighbor hopping element 
with lattice spacing a along the ID ring. Ao = a /2a de- 
scribes the strength of the Rashba SO interaction. <f> is 
the magnetic flux through the ID ring. In our numerical 
calculations, e/ = if disorder potential at each site is 
absent, otherwise, the lattice on site energy e/ is gener- 
ated by a uniform distribution [— w/2,w/2] of disorder 
strength w. Because of the periodic boundary condition, 
c jv+r a ( c N+i,a) is identical as c\ a (ci vCr ). In addition, 
two semi-infinite ID leads are attached to the ID ring to 
make sure that the upper and lower arms have the same 
length. The influence of these two leads is taken into 
account through self-energy terms. The conductance is 
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calculated numerically as outlined in Ref. 

For a certain Fermi energy and Rashba SOI strength, 
the conductance of 3-3 rings array exhibits the periodic 
oscillation as a function of the magnetic flux <I>. FFT 
analysis displays AB, AAS, and other higher frequencies 
as shown in Fig. [2^a). However, as we will show below, 
with the increasing of the disorder strength w, the AB 
oscillation is suppressed and only the time reversal AAS 
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FIG. 3: The AAS and AB oscillation amplitudes versus Ao 
curves are drawn in red and black solid line respectively, for 
a single ID ring (a) under disorder average with w — 0.8, 
Ef = —1.8 and (b) under Fermi energy average. The AAS 
oscillation amplitudes versus Ao at E = — 1.8 (black dotted 
line) has also been plotted for comparisons in (b). Other 
parameters are set as N = 100, to = 1. All the energies are 
measured in the units t . 

Fig. [3£a) shows the AAS oscillation amplitude versus 
Ao curves (red solid lines) and the AB oscillation ampli- 
tude versus Ao curves (black solid lines) under disorder 
average with disorder strength w — 0.8. Compared with 
the AAS oscillation, the AB oscillation has been sup- 
pressed greatly under the disorder average. 

While Fig. [3fb) shows the AAS oscillation amplitude 
versus Ao curves (red solid lines) and the AB oscillation 
amplitude versus Ao curves (black solid lines) under av- 
erage over incident electron energy (Fermi energy). In 
this case, the AB oscillation has not been suppressed. 
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Besides, although these two kinds of averages lead to the 
similar AAS oscillation shapes in the conductance, the 
oscillation amplitude under the disorder average is al- 
most 20 times smaller than that of energy average case. 
The suppression of the amplitude is due to the fact that 
electron will be localized at large disorder strength. 

In these plots, no trace of half oscillation, which ap- 
pears in the experimental results [3], could be found. 
However, from the AAS oscillation amplitudes versus Ao 
for a single ID ring at E = —1.8 [dotted line in Fig. EJb)], 
we notice that the so-called " half oscillation" appears in- 
deed as the experiments 15] . In comparison of the AAS 
oscillation amplitude versus Ao for E = —1.8 and under 
Fermi energy average, it is shown that at the strength 
of SOI where the half oscillation happened, the following 
oscillation amplitude differs a phase shift n with the en- 
ergy averaged oscillation amplitude. When Ao is smaller 
than 0.145 (before half oscillation) , the peak and dip fea- 
tures of these two curves match well with each other. 
But when Ao is larger than 0.145, the peak feature of one 
curve corresponds to the dip feature of another. These 
conclusions can also be inferred from the experiment [TBI] . 
Therefore, we conclude that half oscillation only happens 
at some certain Fermi energy (carrier density). But this 
behavior disappears after ensemble average (disorder or 
Fermi energy average). 

From the FFT spectrum of AAS oscillation amplitude 
versus Ao curves [red curves in Fig. 3(a) and 3(b)], we 
find a frequency peak located at A^ 1 = 31.25. From the 
theoretical time reversal AC frequency a^ 1 — 2rm* /h 2 , 
we get Aq 1 = N/tt. Thus the theoretical time reversal AC 
frequency should be A^ 1 = 31.83 for N=100. Comparing 
these two values, our numerical result agrees well with 
the theoretical value. 

Finally, in order to substantiate our conclusion, we plot 
the conductance versus Ao in absent of magnetic field 
for several cases: without disorder, with disorder aver- 
age and with energy average (Fig. |4]). From their FFT 
spectra, it is seen that different oscillation frequencies are 
included when disorder is absent [Fig. Ufa)]. With strong 
disorder strength, the AB type component is suppressed 
[Fig. 4(b)]. Thus, we are sure that the time reversal AC 
effect have been successfully demonstrated in such a sys- 
tem. However, the case of Fermi energy average is quite 
different from disorder average. In the FFT spectrum of 
energy average [Fig. 0Jc)], other oscillation frequencies 
(two half circle interference ...) can not be smeared out 
under energy average. The AB oscillation amplitude at 
zero magnetic field becomes smaller with increasing dis- 
order strength, but it still remains strong under Fermi 
energy average, which supports our conclusion in Fig. [3] 

In summary, we have numerically demonstrated the 
time reversal AC interference effect in mesoscopic rings. 
In the absent of disorder, various AC interference pat- 
terns are included in the rings. Nevertheless, most inter- 
ference modes are weak against disorder, except the time 
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FIG. 4: The conductance of a single ID ring versus Ao without 
magnetic field, (a), (b) are calculated from Ef = —0.05 with 
disorder strength w = 0, 1.5 respectively, (c) is calculated 
from Fermi energy average without disorder. All the corre- 
sponding FFT spectrum are shown on the right side. The size 
of the ring N=100. The theoretical time reversal (AAS type) 
AC frequency is at A^ 1 = 31.83. In our FFT spectrum, the 
AAS type AC frequency is located at Aq" 1 = 29.30. 



reversal mode. 
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